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For the thermoelectric properties of NaxCo02, we analyze the effect of local Coulomb interaction 
and (disordered) potential differences for Co-sites with adjacent Na-ion or vacancy. The disorder 
potential alone increases the resistivity and reduces the thermopower, while the Coulomb interaction 
alone leads only to minor changes compared to the one-particle picture of the local density approx- 
imation. Only combined, these two terms give rise to a substantial increase of the thermopower: 
the number of (quasi-) electrons around the Fermi level is much more suppressed than that of the 
(quasi-) holes. Hence, there is a particle-hole imbalance acting in the same direction as a similar 
imbalance for the group velocities. Together, this interplay results in a large positive thermopower. 
Introducing a thermoelectric spectral density, we located the energies and momenta regions most 
relevant for the thermopower and changes thereof. 



PACS numbers: 71.27.+a,71.1.Fd 

The efficient conversion between different forms of en- 
ergy represents a primary task in technical applications. 
Thereby the arguably most valuable form is electrical 
energy. On the other hand, heat energy is widely used 
in power plants to produce mechanical energy which is 
then further converted to electricity. The reason for the 
widespread application of this multiple-step process is 
based on the lack of an efficient direct alternative. 
A direct heat-electricity conversion is possible by the 
thermoelectric effect or thermopower. The (dimension- 
less) thermoelectric figure of merit ZT has been estab- 
lished as a benchmark to compare different materials. 
Values of ZT ^ 3 appear to be necessary to trigger an 
extensive industrial use of thermoelectrics^ ^ which could 
then not only replace traditional methods of energy pro- 
duction but could also increase the efficiency of existing 
technical implementations such as cars^ or power plants 
via the conversion of excess heat to electricity. 
The development of novel thermoelectric applications fol- 
lows two routes: First, there is the search for new materi- 
als with a large intrinsic ZT which can either be achieved 
by optimizing features of already known compounds (via 
e.g. bandstructure engineering^) or by the discovery of 
new classes of materials with peculiar properties (e.g. 
Kondo systems^ ^). Second, one tries to diminish phonon 
losses of the thermopower which can be done by avoiding 
collective excitations of the crystal by a specific structural 
arrangement (phonon engineering) such as heterostruc- 
turing^ ^ . 

Opposed to other transport properties, as for example 
the electrical conductivity, the magnitude of the thermo- 
electric effect depends on an asymmetry between the two 
charge carriers holes and electrons. This asymmetry usu- 
ally manifests itself in different velocities and/or spectral 
densities in the vicinity of the chemical potential. 
In this paper, we investigate the thermoelectric prop- 
erties of Naa;Co02 taken as a representative compound 
where the already favorable single-electron features are 



enhanced by the combination of strong electronic corre- 
lations and disorder. We extend the analysis of a previous 
work^ exploring the dependence of the thermopower on 
various parameters, one at a time as well as combinations 
of two of them, e.g. disorder and correlation strength at 
the same time. 

As a non-interacting starting point we used the local den- 
sity approximation (LDA) result from Ref. 10 fitted to 
a tight-binding bandstructure^. Already at the level of 
LDA the thermoelectric performance of Naa;Co02 was 
found to be high and was understood in terms of the 
pudding mold shape of the decisive aig band^. Here, 
we want to analyze the effects of electronic correlation 
in this compound and therefore we have combined LDA 
with dynamical mean field theory (DMFTj^^. The 
LDA+DMFT approach is described in Refs. [TiHTTl and 
we proceed as in Ref. TS' to compute the thermoelectric 
response. We focused on the low-energy excitations which 
give the main contributions to the thermopower. Ear- 
lier LDA+DMFT studies on the spectral properties of 
Na3,Co02 can be found in Refs. IMSSl 
In Fig. [ij we show a summary of the main results de- 
scribed in this study: Considering neither disorder nor 
correlations, the charge carrier properties are described 
by the "pudding mold" band which favors electrons in 
terms of spectral density, but favors holes with respect 
to group velocity. The latter effect is larger so that hole 
transport^^ dominates and we get a positive value of the 
thermopower Figjl] (left panel). Introducing disorder 
to the system increases the scattering rate of both charge 
carriers thus diminishing Figjl] (middle panel). If we 
additionally take into account electronic correlations, the 
hole spectral density is recovered and can be even over- 
compensated whereas the number of electrons is further 
decreased which overall yields a high thermoelectric re- 
sponse, Figjl] (right panel). For a more quantitative pic- 
ture, refer to Fig.[2| where the change of S in the param- 
eter space of correlation U and disorder Ae is visualized. 
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FIG. 1: (Color online) Qualitative effect of disorder and electronic correlations on the charge carrier^^ and the thermopower 
in NaxCo02. On the vertical axis, energy in arbitrary units is shown with the Fermi level Ef and electrons above and holes 
below. The horizontal axis schematically indicates the movement through the material with velocities and charge carrier density 
according to the length and width of the arrow, respectively. The presence of disorder and electronic correlations is visualized 
by scattering centers Ae and the encirclement of the charge carriers accompanied by U, respectively. 
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FIG. 2: (Color online) Change of the thermopower S in the 
parameter space of local Coulomb repulsion U and disorder 
potential Ae for Nao.7Co02 at T = 290 K. The highlighted 
regions indicates realistic parameter values for this material. 



It can be seen that drastic changes happen for example 
for Ae = 0.55 eV where an increasing U leads to an en- 
hancement of S by 200%. Moreover, for even larger values 
of Ae and U = 3.5 eV we found the maximal values of 
S 80 /iV/K in the tested parameter space. Note that the 
system shows insulating tendency for similar values of 
Ae > 0.55 eV and U = 0. This implies that in correlated 
materials disorder effects may be exploited to enhance 
the thermopower. 

The structure of this paper is as follows: In Sec. [Tj we in- 
troduce the Hamiltonian of the model and describe how 
standard LDA+DMET can be expanded by means of co- 
herent potential approximation (CPA) to include disor- 
der. In Sec. [TTj we depict the linear response formalism 
to compute the thermopower as a functional of spec- 
tral density, group velocity and temperature. Both sec- 
tions present the necessary details for following our result 
which were naturally missing in the short paper Ref. |9] 



ally to the lattice structure, the employed tight-binding 
approximation is explained in detail. Confirming previ- 
ous results from Ref. 11, the third nearest neighbor hop- 
ping appears to play a key role for the transport proper- 
ties of this triangular compound. Disentangling the spec- 
tral contributions and the self energy from the two non- 
equivalent lattice sites due to disorder, it is shown that 
one of the two sites exhibits a much larger correlation 
effect than the other. A comparison of low energy spec- 
tral properties to experiment concludes this section. In 
SecjTVj the main results for the thermopower are pre- 
sented. The method of visualization of spectral contri- 
butions to the thermopower introduced in Ref. 9, is de- 
scribed in detail as well as new ways to simplify the theo- 
retical thermoelectric analysis. In an effort to distinguish 
between the parameters temperature, disorder and elec- 
tronic correlation, we investigate their influence on the 
thermopower with the respective other parameters fixed. 
It is revealed that the system is very sensitive to small 
changes of disorder. The impact of correlations on the 
transport properties is directly connected to the strength 
of the disorder which can be seen in a detailed spectral 
analysis. Towards the end of this section, the dependence 
of the thermopower on a change of the sodium doping is 
delineated. Increasing the doping within the investigated 
doping region appears to increase the thermoelectric ef- 
fect which is also supported by experimental results. Fi- 
nally, we give a conclusion in Sec. |Vl 



Sec. Ill provides an insight to the structural and spec- 
tral properties of NaxCo02 missing in Ref. 9, Addition- 



I. LDA-hDMFT APPROACH INCLUDING CPA 



We aim to model NaxCo02 with binary disorder, 
which corresponds to consider two different sites a, 6 with 
an occurrence of Xa^x^ = 1 — x^. In comparison to 
the standard one-band Hubbard mode^, the disorder 



Hamiltonian 



Hdh 



U ^ rii^nn ^ Ae^rii (1) 



ieb 



includes a term contributing only on the site 6 by a dis- 
order potential As. Here, and q^t create and an- 
nihilate an electron on site i with spin a, respectively, 
riia- := cl^Cia- is the electron density on site i for the spin 
cr, U is the local Coulomb interaction. 
For disordered syste ms the Coherent Potential Approx- 
imation (CPApl2Il provides a way to treat inequivalent 
sites in the same spirit as DMFT. In the CPA, the sites 
a, b are assumed to be surrounded by a self-consistently 
determined effective medium. Both sites experience the 
same dynamical mean field, i.e. see the same (averaged) 
environment given by an impurity Greens function Q^^ 
as defined in Ref. 13 or equivalent ly the same (coherent) 
potential represented by a self energy S. Even though 
the mean field Qq^ is the same, the disorder manifests in 
two different on-site energies = 0, £5 = Ae on the sites 
a and 6, respectively (Fig. |3|. 

In DMFT the dynamical mean field Qq^ determines an 
action 



r^rV+(rOV^.(rOV^+.(r)V^.'(r). 



(2) 



where ip are creation and annihilation operators, 
respectively, a is a spin index, {3 is the inverse tem- 
perature and uOrn are Matsubara frequencies (see e.g. 
Ref. nil). The inclusion of (disorder) CPA into DMFT 
is straightforward (see Fig. [s] for a schematic visualiza- 
tion): Assume for the moment that the effective medium 
QQ^{iuJm) is known. Then, the local Greens functions 
Ga{T = {)^iujrn)^ Gb{r = O^iuJrn) cau be determined 
via the action A from ([2| by replacing the standard 
Qo^i^m) by {QQ^{oJm) — l^a,6) with different on-site en- 
ergies Sa = 0, £6 = A^. The local Greens functions Ga, Gb 
are the local interacting propagators of electrons starting 
at a, 6 and returning to the same site a or 6, respectively. 
Then, the total Greens function G{r = O^iuJm) is given 
by the weighted mean 

G{iUm) = XaGa{iuJm) + (1 " Xa)Gb{iuJm) (3) 

such that both propagators G^? G5 are taken into account 
according to the stoichiometric appearance of the corre- 
sponding site Xa or x^ in the lattice. 
For the CPA expansion of the DMFT, the effective dy- 
namic medium QQ^{i(jJm) is determined self-consistently. 
For a given Qq^^ we obtain the local Greens functions 
Ga^Gb by two Hirsch-Fye quantum Monte Carlo (HF- 
QMC) calculations with different on-site potentials = 
0, 65 = As. Then, the weighted G{i(jJm) from Eq. ([3| can 
be used to compute the new self energy T^{uj) or equiva- 
lently the new effective field Qq^ via Dyson's equation 
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FIG. 3: (Color online) Schematic description of the dynami- 
cal mean field method (DMFT) compared to its combination 
with the coherent potential approximation (CPA). Top: In 
standard DMFT the lattice problem is replaced by a site in 
an self-consistently determined dynamical mean field. Bot- 
tom: For the CPA, two sites with different on-site potential 
are connected to the same dynamical mean field. 



With this approach we can include both correlations and 
disorder in a dynamic mean field scheme. 



II. LINEAR RESPONSE AND DMFT 
EXPRESSIONS FOR THE THERMOPOWER 

In linear response, the thermopower is defined as a 
ratio of two correlation functions Ki and Kq 



S 



kBKi 
T Ko' 



(5) 



Here and throughout, T denotes the temperature. The 
determination of the correlation functions usually re- 
quires a representation of current operators in terms of 
the electric field operators The electrical current 

in imaginary time reads 

i\e\ 

j{r,r) = ^ [^^{r,r)\/^{r,r) - \/^^{r,r)^{r,T)] , 
2m ^ 

(6) 

which can be found in many textbooks^^. On the con- 
trary, the representation for the heat current has many 
subtleties and is derived via a continuity equation follow- 
ing Ref. M 



0(r,r) = 



2m 



(7) 



(4) 



where the dot denotes an (imaginary) time derivative. 
Then, the electrical-current-electrical-current correlation 



function evaluates to (as vertex corrections are absent in 
single-band DMFT) 



dm 

duj 



(8) 



with V the unit cell volume, hVe = de/dk the group ve- 
locity of the charge carriers, A(k^uj) the electron spectral 
function and f{uj) the Fermi function. Note that Kq = a 
has units of a conductivity [A m~^]. On the other 
hand, the heat-current-electric-current correlation func- 
tion reads 



duo 



(9) 



where Cgv = 11600 K/eV and Ki has units [A J ^Km 
Combining Kq and Ki leads to the thermopower 
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1^1 Efe |Ve(fc)|2/dw A2(fc^^) 



duj 



(10) 



Note that this result is expressed in the correct units 
[V/K] = [JA-^s-^K] {p times the fraction of and 
remains dimensionless) . 

For the integral in the denominator of Eq. (10), i.e. Kq^ 
the contributions with respect to k and uo enter addi- 
tively, i.e. an asymmetry of \V£{k)\ with respect to k or 
A{k^uj) with respect to (/c,cj) is not important. The inte- 
gral in the numerator of Eq. ([lO| (i^i) has an additional 
uj and since df/duo is negative, \Ve{k)\^ A^{k^uo)ujdf / duo 
is positive for < (holes) and negative for (x; > (elec- 
trons Consequently, asymmetry of the group velocity 
\Ve{k)\ as well as the corresponding asymmetry of the 
spectrum A(k^uj) are decisive for a high thermopowei^^. 
In our model material NaxCo02 this is accomplished 
by a low group velocity V£{k) for electrons, whereas 
holes move much faster through the system. For the 
band structure this means we have a flat region directly 
above > the Fermi level at = and a steep slope 
directly below cj < 0. This is the feature usually denoted 
by "pudding mold" (cf. Ref-E]). 

From Eq. (|8| it becomes clear that Kq is always positive, 
since df / duo is negative for all uo and the other quantities 
only enter quadratically. At the same time Ki may 
be either positive or negative depending whether the 
hole contribution < or electron contribution > 
dominate. This reflects obviously in the sign of the 
thermopower S. 
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FIG. 4: (Color online) Top: Unitcell of NaCo02 according to 
the crystal structure from Ref. 37 We additionally visualize 
the Co and Na bondings as well as one distorted O octahe- 
dron. The actual size of the atoms is different from the hard 
spheres and the distance between the Co and Na layers is 
larger than shown. For the Nao.7Co02 system, we assume 
that only 70% of the Na position are occupied. Bottom: The 
cobalt layers and the distorted oxygen octahedra. 



III. STRUCTURE AND SPECTRAL 
PROPERTIES OF SODIUM COBALTATE 

In 1997, Terasaki and coworkerJ^ identified 
NaCo204 (= Nao.sCoOs) as a promising candidate 
for thermoelectric applications. Other publications con- 
firmed the results of a positive thermopower of the order 
of 100 /iV/K even for other sodium content J^^^^^. With 
appearance of superconductivity in the H2O intercalated 
compounds^^ the material NaxCo02 once again 
proved to exhibit surprising features. 
Recently, Kuroki et al}^ argued in favor of band 
structure properties playing the most important role 
in determining the thermopower in NaxCo02. In their 
study, Kuroki and co-workers used Boltzmann's equa- 
tion to compute the Seebeck coefficient S{T). Here, we 
go beyond that approach and propose an "ab initio" 
method to compute 5* by means of linear response 
theory including correlations by DMFT and disorder by 
CPA (cf. Sec, [m and [H respectively). As representative 
compound, we choose a doping of XNa = 0.7, i.e. 
Nao.7Co02. 



Lattice structure 



In Fig. |4j we plot the crystal structure of the system. 
Sodium cobaltate NaCo02 consists of alternating lay- 
ers of sodium and cobaltate C0O2. The Co atoms are 
arranged in a trigonal lattice structure with a nearest- 
neighbor distance of 2.84 A, see Ref.[37l Each cobalt atom 
is surrounded by 6 oxygen atoms in form of a tilted oc- 
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FIG. 5: (Color online) The breakup of the degeneracy of the d 
orbitals: First, the 5 d basis energetically split up in 3 t2g or- 
bitals and 2 Cg orbitals in an octahedral configuration. Second, 
the 3 t2g orbitals split up in one aig and 2 e^/ due to the dis- 
tortion of the oxygen octahedra. The energy differences shown 
in this plot are symbolical. We also visualize the occupancy of 
the orbitals in the different situations, where oxidation state 
of the cobalt atom in the compound is Co^^ (Co^^) if a Na 
partner is (not) present. This additional electron whose pres- 
ence depends on a sodium donor is shown as an arrow with 
different color. 



tahedron. The Na atoms are located in their own layer 
half way between the C0O2 layers. 

The exact lateral position of the sodium atoms within 
their layer and mutual shifts of the C0O2 layers, de- 
pend non-trivially on doping, temperature, etc. For ex- 
ample, Huang^^ reports a change of the Na positions for 
^Na = 0-75 around T = 320 K. On the other hand, 
an earlier work of Yakabe^^ finds three typical phases 
of NaxCo02 as x-^a changes from 0.55 to 1. Since almost 
all of the transport takes place in the C0O2 layers, mu- 
tual shifts of the cobaltate layers only weakly influence 
the thermopower results. We therefore assume the C0O2 
layers to be all equivalent, i.e. there is no structural shift. 
The lateral position of Na in the compound, on the other 
hand, is a delicate topic. Many Na coordinates have been 
reported, mainly the position directly above the Co atom 
{x = O^y = 0) and a positio n whi ch is aligned to the 
O atoms {x = 2/3, y = 1/3^^. Moreover, the for- 
mer reference also suggests movement of the sodium ions. 
For lower temperatures Zandbergen^^ found two different 
types of Co atoms, namely Co^+ and Co^+. We therefore 
assume that the cobaltate system shows binary disor- 
der, even if the sodium atoms are distributed on mul- 
tiple sites. As we are interested in relatively high tem- 
peratures for thermoelectric applications, a possible or- 
dering of the Na ions does not appear to be relevantPSl. 
For simplicity, we assume the Na to be located on the 
(x = 0, 7/ = 0, 2; = 1/2) site directly above the cobalt. 
Thus, we can still assign a Na atom to one specific Co site. 
This corresponds to the same "virtual crystal approxi- 
mation" which has been used for the underlying band 
structure calculatiorP^. 

For a doping of XNa = 0-7, 70% of the sodium sites are ac- 



tually occupied and we consequently have 30% vacancies 
in the sodium layers. The hereby inequivalent Co sites 
in the cobaltate layers are denoted by CoNa (= Co+^) 
and Co Vac (= Co+^) depending on whether they have a 
close sodium atom or not. As we will show below, elec- 
trons of one band with little kz dispersion are responsible 
for transport. We denote the filling of this band at the 
two different sites by UNa and nvac, respectively. Due to 
the disorder, these electrons tend to be at cobalt sites 
CoNa with a sodium partner, by means of simple elec- 
trostatic attraction to the Na+ ionic cores. The strength 
of this tendency is governed by the disorder potential 
As, i.e. the difference of the on-site potential of the two 
non-equivalent sites. Additionally, many-body effects can 
influence the movement of electrons and holes in the dis- 
ordered crystal. 

To determine the relevant energy bands in this material, 
it is vital to understand the crystal field splitting of the 
NaCo02 lattice. In this light, as part of a thought ex- 
periment, the gradual transition between the Co atomic 
orbitals and their bulk material analogs can be useful: 
The 5 atomic d-orbitals can be classified according to 
their radial and spatial dependence of the corresponding 
wave functions. In a Cartesian coordinate system one has 
three planar orbitals \xy), \yz), \zx) denoted by t2g and 
two orbitals |3z^— r^), \x'^ — y'^) denoted by Cg. Refer 
to Fig. [5] for a visualization of these states. The electron 
configuration of single- atom cobalt is [Ar]3d''4s^. Accord- 
ing to Hund's rule, the d-electrons try to occupy all the 
orbitals and maximize the total spin (in this case S=l), 
if the d-orbitals are degenerate. 

Atomic orbitals hybridize if the atom is placed in bulk 
material. However, one can use the picture of atomic or- 
bitals as a starting point to gradually include the influ- 
ence of the crystal structure. In C0O2, the Co atom is 
surrounded by an octahedron of 6 O atoms. The Cg or- 
bitals of Co have a large part of their spatial probability 
density close to these ligands. Since the O ion is nega- 
tively charged, the Cg orbitals pointing towards the oxy- 
gen become less favorable for electrons to occupy than in 
the atomic case. 

Moreover, the O octahedra which surround the Co atom 
are distorted, cf. Fig. |4] One can show, that this does not 
lift the degeneracy of the e^, but, on the other hand, the 
t2g states split up in 1 aig and 2 energetically lower e^/ 
orbitals. These aig and e^/ states do not correspond any 
longer to single t2g orbitals, but are rather linear combi- 
nations thereof. These considerations lead to the orbital 
structure visualized in Fig. [5] 

When we go beyond a single Co site, the picture of local 
orbitals is not fully valid any more, but is rather replaced 
by bands representing a certain dispersion relation. Still, 
we can look for bands in the band structure which have 
predominant aig, e^/, or Cg orbital character and which 
arrange in the same energetic order as the orbitals in 
Fig. [5] 

The band structure from Singfl^ provided the starting 
point for the tight-binding approximation. To mimic dis- 



FIG. 6: (Color online) Tight-binding approximation (solid 
line) and group velocities (dashed lines) for Nao.7Co02. 



order at the level of LDA, Singh used a NaCo204 super- 
cell with twice as many bands due to the doubling of the 
unit cell. Following a different approach, we extract an 
approximation for one single ai^ band of NaCo02 from 
the LDA result and include the disorder in a second 
step at the level of the DMFT algorithm, cf. Sec. [T[ 
Corresponding to atomic valence orbitals, one also finds 
valence bands within a solid state compound. Analo- 
gously, very tightly bound "core bands" and, in the 
ground state unoccupied, "conductance bands" may be 
identified. In cobaltate, the valence bands are the 3 oxy- 
gen p-orbitals and the 3 cobalt d-t2^-orbitals. Due to high 
electronegativity, the p-bands of O are well below the Co 
bands. On the other hand, the sodium s-orbitals belong 
to the conductance bands and are even above the unoc- 
cupied Co e^-bands. The number of bands considered in 
the LDA+DMFT model strongly influences the compu- 
tational effort necessary to solve it. Whereas the decrease 
to the complexity for the tight-binding approximation is 
not decisive, the effort for the DMFT calculations heav- 
ily depends on the number of bands. We aim therefore 
to minimize the number of orbitals, without losing (too 
much) physically relevant information. 
The oxidation state of Co in the compound is +3 in 
NaCo02 and +4 in C0O2, respectively. Thus, of the 9 
valence electrons 6 and 5, respectively, actually fill the 
cobalt d bands. For transport properties a rather small 
energy interval of the order of 0(/cbT) around the Fermi 
edge is important, since electrons in this regime can 
be thermally activated. Only the t2g bands are located 
within that interval (under normal conditions) . In the 
following, we will assume that the hybridization of e^/ 
and aig bands is negligible. With this assumption, the 
crystal field splitting visualized in Fig. [5] guarantees that 
the orbitals remain filled and that, for transport cal- 
culations, restricting to the partly occupied aig orbital is 
sufficient. 

To extract the two effective aig bands from the band 



FIG. 7: (Color online) Density of states N(uj) considering an 
increasing number of neighbors in the tight-binding approx- 
imation using the parameters from Table ^ The dotted line 
corresponds to the DOS computed only with nearest-neighbor 
hopping, whereas the solid curve is computed with all 9 pa- 
rameters. A decisive hopping appears to be the one to the 
third nearest-neighbor, where a significant change of the char- 
acteristics of the DOS can be observed. 



structure of NaCo2 04, Kuroki and Usui made a tight- 
binding fit^, which is shown in Fig. [gJ Since we aim to 
include only one ai^ band in the DMFT, we have to find 
the corresponding single tight-binding band for NaCo02 . 
Though NaCo02 is similar to NaCo2 04, there are sub- 
tle differences: Since the super cell for NaCo204 contains 
two Co atoms the corresponding Brillouin zone has half 
the size of the one for NaCo02 in z direction. Thus, 
the single aig band will not match exactly one of the 
two bands for NaCo204 of the original fit. The lower 
band of the two highlighted bands mainly belongs to the 
cobalt site in NaCo204 with a sodium partner. As an ap- 
proximation, we therefore take the tight-binding param- 
eters for this band, but use the (larger) Brillouin zone 
of NaCo02 in DMFT. The final aig band is shown in 
Fig. |6] together with the corresponding group velocities 

In the following, we discuss the applied tight-binding ap- 
proximation: Let the aig band have the dispersion re- 
lation e^^^{k) given by the data points which are the 
output of the LDA calculation. We now aim to replace 
e^^^{k) by an analytic expression, which then can be eas- 
ily used to obtain the group velocity Ve^^^{k). To this 
end, we make a tight-binding ansatz with 9 parameters 

9 #NN 
n=l j=l 

where #NN denotes the number of nearest-neighbors 
of nth order and rj^ enumerates the corresponding nth 
nearest-neighbor position. Due to inversion symmetry the 
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exponential function in (11) can be simplified as 



,LDA 



(k) 



(fc) 



9 #A^A^ 
n=l j=l 



(12) 



Here, ^NN means that the pairs of atoms connected by 
inversion symmetry are counted only once. The first 4 pa- 
rameters ti, . . . , ^4 describe the hopping between nearest- 
neighbors within one Co-plane. The fifth parameter t^o 
measures the probability of Co-interplane hopping from 
a cobalt atom to another cobalt directly above (below). 
The final 4 parameters t^i, . . . ,^^4 describe processes of 
combined interplane and nearest-neighbor hopping. 
The fitted set of parameters^ are listed in Table |lj The 
final tight-binding fit of the ai^ band yields a bandwidth 
of 1.5 eV and is, as expected, predominately dispersive 
in the x-y-plane (corresponding to the cobaltate layers). 
With the analytic dispersion £™(A:) from Eq. (12), we 
may investigate which tight-binding parameters are nec- 
essary for the characteristics of the non-interacting LDA 
density of states 



N{uj) := ^ 

TT 



E 



1 



u ■ 



{k) + iS 



(13) 



with ^ ^ 0+. The functions N{u) considering hopping 
up to 1, . . . , 9 nearest-neighbors are plotted in Fig. [t] (the 
final solid line corresponds to the LDA-DOS). Consistent 
with the results from Kuroki^H we observe that hopping 
up to the third nearest-neighbor is necessary to obtain 
the important upper van-Hove-like peak, which will play 
a role for the explanation of the large thermopower below. 
If we also consider interplane hopping (i.e. tz^ 7^ 0) the 
peak above = is less pronounced. Leaving the 2d 
regime this is the standard behavior for van-Hove peaks, 
which are usually damped in 3d with respect to lower 
dimensions. The final DOS is dominated by two major 
peaks which are approximately 1 eV apart. 
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intra / interlayer 


NN 


values (par. fit) 


ti 


intralayer 


1st 


0.1800 




intralayer 


2nd 


-0.0388 


U 


intralayer 


3rd 


-0.0270 


u 


intralayer 


4th 


0.0004 


tz,0 


interlayer 




-0.0180 


tz,l 


interlayer 


1st 


-0.0049 


tz,2 


interlayer 


2st 


-0.0016 


tz,3 


interlayer 


3rd 


0.0011 




interlayer 


4th 


0.0005 



TABLE I: Tight-binding parameters obtained by a parameter 
fit to LDA data. 
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FIG. 8: (color online) Top: k-integrated spectra A(l)j) obtained 
by DMFT and maximum entropy method compared to the 
non-interacting DOS N{uj). Bottom: he self energy ^^(icj) 
on the imaginary axis. Note that sodium and vacancy contri- 
butions are resolved. 



B. LDA-I-DMFT self energy and spectrum of 
Nao.7Co02 

With the tight-binding dispersion relation e^^, we now 
run the DMFT algorithm depicted in Sec. [T[ For given 
temperature T, correlation /7, doping XNa, and disor- 
der potential As, we obtain the self energy T^{i(Jrn) at 
the Matsubara frequencies oorn = (2m + l)7r//3, and the 
Greens function G{r). Here, and throughout this work, 
P denotes the inverse temperature in eV~^ defined by 
13 := 11600/T eV~^. From the imaginary part of the 
self energy ^S(ia;) we can extract information about the 
strength of correlation and disorder effects in the system. 
For example, a divergent 



lim ^Ti(iuj) = ±00 



(14) 



indicates that the charge carriers are localized which 
leads to an insulating phase. On the other hand, a corre- 
lated metal usually yields 



lim ^S(zco') = 0, 



(15) 



with a quasiparticle weight Z ^ 1. 

In Fig. [s] (bottom), we plot the self energy over the 
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uieV) 

FIG. 9: (Color online) The two contributions A]sia{(jo) and 
AYac{(^) to the overall spectrum Atot{(^) without the stoi- 
chiometric weighting factor. The filling at the sites CoNa is 
nNa= 0.91, therefore correlation effects are negligible. For 
Covac with nvac= 0.71 on the other hand, the formation of a 
quasiparticle peak at the Fermi edge cj = 0, a lower Hubbard 
band at u ~ —U = —3.5 eV and a upper Hubbard band at 
Lj ^ 0.5 eY can be observed. The center of the two spectra are 
approximately shifted by the disorder potential Ae = 0.55 eV 
with respect to each other. 



Matsubara frequencies ujm for our set of initial parame- 
ters 

{T = 290 K,U = 3.5 eV, XNa = 0.7, Ae = 0.55 eV}. 

(16) 

Note that here we employ similar parameters as in a pre- 
vious LDA+DMFT calculation for Nao.7Co02^^ which 
did, however, not analyze transport properties'^. At the 
first glance, the overall ^Utot shows characteristics of a 
bad metaP^, since in the whole temperature range con- 
sidered it neither diverges nor converges to for ^ 0. 
Hence, the separation of correlation and disorder effects 
on needs further investigation. As a first step, we 
can also visualize the contributions of the two sites CoNa 
and Covac to the self energy ^Stot- Following Sec. [T| 
these contributions are computed via Dyson's equation 
with the corresponding Greens function, i.e. 

SNa= Go'^^e-G^l (17) 

Svac = Go^ - Gyl^, (18) 

Stot= go'-G-\ (19) 

where G = ^Na^Na + ^VacG^Vac- We observe that the self 
energy ^I^Na for the sites CoNa is metallic, the other con- 
tribution ^Svac for CoNa showiug insulating tendency. 
Aquivalent to the self energy ^I](za;) connected to the in- 
fluence of correlation and disorder, the Greens function 
G(r) provides information about the many-body excita- 
tions of the system, as the k-integrated spectrum A{uj) is 
given implicitly by 
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FIG. 10: (Color online) Spectral function A{k,uj) for T = 290 
K obtained by a Taylor fit of the DMFT self energy. The 
color code of the values of A{k, uo) is given by the colorbar 
on the right-hand side. The dashed line corresponds to the 
non-interacting tight-binding dispersion relation. 



/ r(/i-u;) 

To calculate the spectral functions, we employ the max- 
imum entropy method (see for example Ref. I43|44p . In 
Fig. [s] (top), we plot the total spectrum Atot{^) for the 
diig band and compare it to the contribution of the two 
sites CoNa,Covac as well as to the non-interacting density 
of states (DOS) N{uj). Note that we do not show here the 
lower Hubbard band at ~ —U = — 3.5eV. At the first 
glance, the changes between the N{uj) to Atot(^) are not 
significant. But since two contributions ^Na, ^vac with 
two different electron fillings 

nNa = 2 • 0.91, nvac = 2 • 0.71 (ntot = 2 • 0.85 = 1 + XNa) 

(20) 

mix, a separate investigation of the two spectra is use- 
ful (note that nNa^Na + ^Vac^Vac = ^tot). On CoNa sites, 
the diig band is almost filled, whereas on Covac the value 
is closer to half filling, where correlation effects are ex- 
pected to be stronger. In fact, as can be seen in Fig. the 
formation of a quasiparticle peak at = and a lower 
Hubbard band at ^ —U = — 3.5 eV in the spectrum 
Avac for the sites Covac can be observed. Especially the 
Fermi liquid behavior for small cj, may be a decisive ingre- 
dient for higher thermopower, since the spectral weight 
in the thermally activated energy interval is crucial as 
already discussed in Sec. [llj This analysis indicates that 
disorder enhances the effect of correlations in the com- 
pound driving the sites Co with Na vacancy, closer to 
half filling. As a consequence, the quasiparticle peak for 
a given correlation U is more pronounced than in the case 
without disorder. 

Analogously to the filling, one can also investigate the 
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approximate renormalization of the two bands 

ZNa = 0.88, Zvac = 0.36 (Ztot = 0.48), (21) 

which was obtained from ^Ti{iujm) by using the tan- 
gent of the two Matsubara frequencies closest to the 
origin. This values should be considered with caution, 
since the disorder perturbs the Fermi liquid behavior at 
both sites CoNa and Covac- However, as Fig. [9] indicates, 
Nao.7Co02 can be interpreted as a system consisting of 
sites with large correlations Covac, and sites CoNa where 
correlation effects are weak. 

Apart from the k-integrated spectrum A{uj)^ we are also 
interested in the k-resolved quantity A{k^uj)^ that is in- 
serted in Kubo's formula (10). Therefore, we need an an- 
alytic continuation T^{iuj) S(cj), which is, in general, 
non-trivial. The behavior of S(a;) further away from the 
Fermi level is not important for transport properties and 
a Taylor fit is hence sufficient. With the self energy T^{uj) 
the spectrum is calculated as 



A{k,uj) 



1 



(22) 



where /i is the chemical potential from the DMFT calcu- 
lation. 

A visualization of the approximation for A(k^uj) com- 
pared to the tight-binding band can be found in Fig. [lO] 
Let us remark here that it cannot be expected that 
the electronic structure or spectra match the experimen- 
tal results in a large energy interval around = for 
a model without the egf bands. However, ARPES data 
from Yang^^ indicates that the e^/ bands, which appear 
to cross the Fermi edge, actually do not contribute to the 
Fermi surface, i.e. are shifted down in energy by many- 
body effects. Thus, a comparison to experimental data 
may well be accurate in the vicinity of the Fermi level. 



In Fig. 11 , we show the numerically obtained spectra and 
the corresponding ARPES data^^. For a direct compari- 
son, all numerical results are already multiplied with the 
Fermi function /(cj, /3) = (e^^+i)-^ with T = 290K as in 
the experiment (we thus assume that the spectra do not 
change qualitatively as Tnum = 290 Texp = 40 K). 
On the left-hand side the k-resolved spectrum A{k^uj) is 
plotted for various values of k along the /c-path T ^ K.ln 
the inset, we give the corresponding experimental results. 
The reason for the good agreement is the large diig char- 
acter of the experimental spectra. On the right-hand side 
of Fig. 11, we compare the k-integrated spectrum A{uj). 
To simulate the e^/ background we add a linear function 
for uj < 0.2 eV. Again the data coincide qualitatively, as 
the inset with the experimental results shows. 
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FIG. 11: (color online) Top: Comparison of the spectrum 
A(k,uj) to ARPES results, inset taken from Ref. 45). In ex- 
periment, the peaks correspond to aig excitations. Bottom: k- 
integrated spectrum A(uj) and corresponding experimental re- 
sults (inset taken from Ref. I45p. The two slopes in the inset 
are spectra for different photon energies. For comparison, a 
linear background signal for cj < —0.2 eV has been added by 
hand to the numerical results. 



IV. THERMOPOWER 



In order to understand the influence of temperature, 
disorder, correlations and doping on the self energy and 
the thermopower, in this section we will investigate their 
effect separately with the other parameters kept fixed. 
For disorder and correlations, however, it is revealed that 
their effect strongly depends on the respective other pa- 
rameter. We therefore discuss combined effects of disor- 
der and correlations in an own subsection. Before start- 
ing with the first parameter, let us introduce some helpful 
quantities which will be helpful in the analysis. In Sec.[ll| 
we derived the expression S{T) ~ Ki/Kq for the ther- 
mopower. The pivotal quantity which we now focus on, is 
Ki depicted in Eq. ([9| and therein especially the kernel 
as a function of (/c, u). In most of the cases, changes of Ki 
directly reflect on the thermopower S. If we separate the 
factors in the kernel for the initial parameters (16), we 



arrive at contributions visualized in Fig. 12 Investigating 
these contributions can lead to a deeper understanding 
of the important processes in thermoelectric transport. 
In order to factorize the total kernel /Ctot, we define 5 
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auxiliary functions of k and w 

df{uj,T) 



/Ci(/c,cj) := 



duj 



df{oj,T) 
K.2{k,uj) := — — — w, 

K:s{k,uj) := \Ve{k)\\ 

lC,{k,oj) := \Ve{k)\'^^^^ cv 

/C5(fc,w) :=A^{k,uj)^^^^u;, 



(23) 

(24) 
(25) 
(26) 

(27) 



TKo V 



(28) 



where we included all prefactors in /Ctot such that S = 
/ dujJCtot{k,uj). The total kernel /Ctot may thus be 
interpreted as a "thermopower spectral density" deter- 
mining the contribution to 5* in (/c, cj)-space. In the fol- 
lowing we restrict ourselves to in-plane transport, i.e. we 
set I Vef^^ =v'^ = vl+ vl. 



In Fig. 
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a visualization of /Ci, . . . , /C5, /Ctot for {T = 
3.5 eV, = 0.55 eV} is shown. The deriva- 



290 K, U 

tive of the Fermi function /Ci constitutes an energy inter- 
val around the Fermi edge, where electrons are thermally 
activated. There will be no electronic transport outside 
of this interval. The different sign of the charge carri- 
ers electrons (negative) and holes (positive) is accounted 
for by the factor in /C2. Thus, Fig. [l2|top, right) shows 
the charge carriers activated for transport with respect to 
sign. On the other hand, the group velocity /C3 is a quan- 
tity which only depends on the band structure. As a con- 
sequence, /C3 will be large whenever the slope of \/e{k) 
(middle, left). Combining the latter 

(mid- 



12 



12 



is steep, cf. Fig. 

three contributions leads to /C4 shown in Fig. 
die, right). We observe that the values are still distributed 
symmetrically with respect to the Fermi level. The effect 
of the temperature and spectrum without the group ve- 
locity is encoded in /C5. There appear to be more contri- 
butions from the electrons uj > (Fig. [12] (bottom, left)). 
Finally, the thermopower spectral density /Ctot indicates 
the regions of {k^uj) space which contribute to the ther- 
mopower, cf. Fig. [12] (bottom, right). The function /Ctot 
yields both positivefholes) and negative (electrons) con- 
tributions to the thermopower, which may in principle 
also annihilate each other. In fact, that is the usual be- 
havior for most materials where the thermopower S is 
therefore only of the order of ±0(l)/iV/K. In our exam- 
ple, the electron contributions are significantly smaller 
than hole contributions, and this inbalance reflects itself 
on a positive thermopower S. 



A. Effects of the temperature 

Before considering the numerical results of the model 
system Nao.7Co02, we discuss the effects of temperature 



to the thermopower 5* by a rather general approach ana- 
lyzing the thermal dependence of the kernel /C5(/c, u). Let 
us assume for the moment the extreme case that the spec- 
trum A{k^uj) ~ A{k)Q{—uj)^ i.e. only holes contribute to 
the thermopower. Then, we can solve the integral 



J — cx 



duo 



UJ 



(29) 



to obtain an approximation of the thermal dependence 
of the thermopower S: 



m = f{oj,/3)uj\ 

ln(2) 



dujf{uj,/3) 



11600 



T [eV]. 



Analogous considerations for the conductivity a = Kq 
yield a factor P which cancels with P from the definition 
of 6* = k^Pe'^KiK^^ . Thus, a linear dependence S (xT 
represents a first approximation for the thermopower. 
Since extreme asymmetries as assumed in Eq. (29) will 



typically not persist for larger T the linear dependence 
will be reduced for higher temperatures. 
For low temperatures T, the spectrum A inside of the 
thermally activated energy interval will become symmet- 
ric with respect to cj, i.e. A{k,uj) ~ A{k). Then, the 
integral 



HP) 



f 

J — c 



duj 



uj = 



vanishes, since the integrand is an odd function. Thus, 
for small T ^ 0, the thermopower S is expected to con- 
verge to 0. 

After these general considerations, we now return to 
Nao.7Co02: In Fig. [13) we plot the imaginary part of the 
self energies ^T^{iuj) on the Matsubara axis and the k- 
integrated spectra A{uj) for various temperatures T. The 
high energy tail of ^E(zco') is determined by a sum ruld^^, 
and is therefore almost temperature independent. In the 
vicinity of = 0, we observe a variation of ^S(zco') with 
respect to temperature. This has several reasons: First, 
the highest Matsubara Um moves closer to a; = by def- 
inition and therefore we obtain a finer resolution. On the 
other hand, incoherent effects due to electron-electron 
scattering decrease with lower temperature. 
The k-integrated spectrum plotted in Fig.[l3]^top), shows 
less spectral weight of A{uj) in the vicinity of the Fermi 
edge for increasing temperature T. This is a usual effect 
of higher temperatures: the dampening of narrow peaks 
in the spectrum. In addition, there are no indications for 
a phase change to an insulator of the material in the in- 
vestigated temperature interval 145-1160 K originated in 
the electronic structure, since the qualitative spectral dis- 
tribution remains more or less unchanged. However, the 
result for /3 = 20 [eV~^] is puzzling, because it shows a de- 
viation of the general temperature behavior with respect 
to the other results. The same behavior is also observed 
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FIG. 12: (color online) Contributions /Ci, . . . ,/C5,/Ctot to the thermopower S for T=290 K according to Equations (23)-(28). 
Positive values are emphasized by bright (yellow) domains, negative contributions are visualized by dark (red) domains. The 
units of the kernels /Ci , . . . , /C5 are arbitrary, but domains with brighter or darker scale (intenser red or blue colors) correspond 
to larger positive or negative values, respectively. As a dashed line, we also plot the LDA band structure. The solid line shows 
the isovalues ±10^ /SVK~^eV~^ for the thermopower spectral density. 
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FIG. 13: (color online) Top: k-integrated spectra A(uj) for var- 
ious temperatures in eV~^. The spectral weight around the 
Fermi edge is decreasing with respect to temperature T. There 
is no major change of ^S(icj) with respect to temperature, 
except for Matsubara frequencies uom close to the cj = 0. Bot- 
tom: Self energy ^^(icj) on the Matsubara axis for various 
temperatures. 



for other doping x-^a = 0.6, 0.8. Most probably this re- 
sult is connected to the numerical error in the employed 
Maximum Entropy Method (MEM) as two spectral peaks 
merge into one. 

In Fig. [T4j we show how the k-resolved spectra and the 
thermopower spectral density /Ctot changes with temper- 
ature. We observe that the spectral weight is in fact de- 
creasing in the thermally activated energy interval. Thus, 
the linear approximation for S{T) has to be adapted for 
the spectra of Nao.7Co02. The linear increase will be 
dampened with increasing temperature leading to a flat- 
tening of S{T). This can also be seen in the final results in 
Fig. 25 , where we plot the thermopower S{T) for various 



B. Effects of the disorder 

Separating the effects of correlation and disorder is a 
highly non-trivial task. Generally, correlations seem to 
have a larger impact on the thermopower S than disor- 
der in Nao.7Co02. With the aim of understanding, we 
first set U = eV and gradually increase the disorder 
potential As = 0^1 eV. In the next section the oppo- 
site case, namely no disorder and increasing correlation, 
will be investigated. 

The probability for a charge carrier to be scattered from 
one site CoNa^Covac to the other is obviously propor- 
tional to the doping XNa- On the other hand, the strength 
of the scattering event is determined by the disorder po- 
tential As. In fact, Ae can be interpreted as an "electron 
affinity", since a larger Ae will lead to more electrons 
occupying the sites CoNa by means of the electrostatic 
attraction of the sodium ions. 



In Sec. IIIB, the importance of the self energy Ti{iuj) on 
the imaginary axis has been discussed. With no correla- 
tion U = 0, we can observe the pure scattering effects 
of the disorder. In Fig. [Tsj the imaginary part of the 
self energy ^T^{iuj) is shown. The value lima;^o ^S(za;) 
gradually increases as Ae = 0^1 eV, which can be 
understood as follows: With the Greens function G := 
G(0) and the Weiss field Qq^ := Qq^{0)^ we make a 
CPA ansatz) 



G 



XNa 



5o 



Ae 



1 - ^Na 



(30) 



This means that the total (local) Greens function consists 
of contributions corresponding to the sub lattice CoNa 
and Covac, where the former is energetically lowered by 
the disorder potential A^. As in the atomic limit, the 
task is to extract the self energy S = S(0) by expressing 
G of Eq. ( 30 ) with the corresponding Dyson equation 



G = 



After some algebra, we find 



(31) 



S = XNa Ae 



(1 - XNa) Ae JR^o + i{l - a;Na) Ae 9C?o 



[l-(l-a;Na) Ae3?go]2 + [(l 



XNa) Ae 95o]^ 

(32) 



To obtain an approximation of the Weiss field ^q^^ , we use 
the Fourier transform of the non-interacting local Greens 
function Go{uj = 0) 



Go 



■I 



V / de 



N{e) 



i7r7V(0), 



(33) 



doping XNa 



where we used the Sokhatsky-Weierstrass theorem and 
N{e) denotes the non-interacting DOS. Note that we 
used the advanced version of the theorem, which means 
we will restrict to the part ujm < of the self energy 
Ti{iujrn)- Now, we assume a constant and symmetric DOS 
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FIG. 14: (color online) k-resolved spectra A(k,u) (left) and corresponding thermopower spectral densities /Ctot (right) for 
various temperatures T. The color code is given by the colorbar. The dashed line on the spectral images is the tight-binding 
fit. The spectral weight around the Fermi edge is decreasing with temperature. Thus, the increase of the contribution to the 
thermopower gets dampened and is therefore not linear. 
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FIG. 15: (color online) The imaginary part of the self energy 
^E(icj) on the Matsubara axis for T = 290 K and t/ = eV 
for various values of the disorder potential As = 0^1 eV. 
A gradual increase of (iu) with respect to As is observed, 
(inset) The values lima;^o ^S(icj) obtained by the approxi- 
mate formula (34) and by DMFT. For small disorder As the 



N{e) ^ N{0)e{e - eo)e{eo - e) with N{0) ^ 0.8 eY'^ 
from F ig. [7| To that end, the principle value term drops 
in Eq. (|33|) and together with Eq. (32) we arrive at 



lim ^S(zco') 



l + [^(l-XNa)A£7V(0)]2 



(34) 



The values of lim^^^o ^S(zco') as a function of Ae are 
plotted as an inset in Fig. [15) and compared to the 
self-consistently determined DMFT+CPA results. The 
DMFT+CPA data indicate quadratic dependence with 
respect to Ae, which is properly approximated by our 
simple considerations. 



approximation is in good agreement with the DMFT data. 



The approximate formula ( 34 ) indicates a direct connec- 
tion between the disorder {^Na, and the self energy 
at low frequency lim^^^o ^S(icj). In fact, the value ^11(0) 
can be interpreted as an inverse time r^g^, the average 
time which an electron spends at a site CoNa oi" Covac^ 
before scattering to the respective other site. 
The increased scattering rate due to disorder is affect- 
ing both Kq and Ki. Numerical results indicate that Ki 
is stronger suppressed than the conductivity Kq which 
corresponds to an overall drastic decrease of the ther- 



mopower S which is visualized in Fig. 16 
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FIG. 16: (color online) The thermopower as a function of 
the disorder potential As without correlation {U — 0) for 
T — 290 K. The inset shows the corresponding values for the 
resistivity p. The magnitude of 5* ^ Ki/Kq becomes smaller 
with larger As, which can be traced back to a decrease of the 
current-heat correlation function Ki dominating the increas- 
ing p = 1/Kq. 
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FIG. 17: (color online) The imaginary part of the self energy 
^Yi(iuj) on the Matsubara axis for different values of the cor- 
relation [/ = ^ 5eV, and T = 290 K (no disorder). The 
actual values of IJ are according to the colorbar on the right- 
hand side, (inset) The values of lima;^o ^^(icj) as a function 
of U. 



C. Effects of electronic correlation 

After considering disorder without correlation in the 
previous section, we now focus on the effects of electron- 
electron correlation and assume that there is no disorder. 
At half filling n/2 = l/2,itis well known that for increas- 
ing correlation [/, a Mott metal-insulator transition can 
be observed (see e.g. Ref. »13;) . 

Away from half filling n/2 < 1/2, an electron occupying a 
certain site simply hops to a neighboring empty site with- 
out the necessity of a double occupancy (which would 
cost the energy [/). Thus, we additionally obtain stan- 
dard conduction channels. An equivalent picture is also 
applicable for holes for n/2 > 1/2 as in our model mate- 
rial Nao.7Co02. Without disorder, we have a mean filling 
of n/2 = 0.85, which is rather far away from n/2 = 1/2. 
For U — 3.5 eV, a renormalization of the non-interacting 
energy band of Z = 0.77 is obtained. Consequently, with- 
out disorder, electronic correlations are only intermedi- 
ately strong. However, it pays to investigate the effect of 
correlation on the self energy to compare with the case 
of combined disorder and correlation. In Fig. [T7| we plot 
the imaginary part of the self energy Zl(ici;) for differ- 
ent values of the correlation U . Several features can be 
observed: First, the maximum max{^S(ico')} ~ U . Sec- 
ond, the limit lim^^o ^S(icj) / and the value depends 
on U . This behavior can be explained by means of Lan- 
dau's Fermi Liquid Theory, which predicts a dependence 
lim^^^o ^^l^*^) a{U)T'^. A visualization of the values 
lim^^o ^^{i^) as a function of U can be found in the 
inset in Fig. [T7[ 

In Fig. [T8| we show the thermopower and the resistivity 
as a function of U with Ae = 0. Larger U appears to in- 




FIG. 18: (color online) The thermopower 5* as a function of 
the correlation U without disorder (Ae = 0) for T = 290 K. 
The inset shows the corresponding values for the resistivity p. 



crease p which reflects on the thermopower S. However, 
the change of S is much smaller than in the case of non- 
zero disorder Ae > below. 

Correlation effects are small without the simultaneous 
consideration of disorder. When we include disorder, one 
of the two sites, namely Covac, will be driven closer to 
half filling. Then, the formation of a quasiparticle peak 
can be observed, which eventually leads to a higher ther- 
mopower. 



D. Effects of combined disorder and correlation 

The pure effects of disorder and correlation on self en- 
ergy and spectra were discussed in the previous two sec- 
tions. Whether the observed characteristics of disorder 
and correlation remain qualitatively unchanged in the 
presence of the other parameter, i.e. correlation or disor- 
der, respectively, is yet to be investigated. We therefore 
fix either U = 3.5 eV or Ae = 0.55 eV and vary the re- 
spective other parameter (T = 290 K and = 0.7). 
In Fig. 19, the imaginary part of the self energy ^I](za;) 
for different values of the disorder potential Ae = ^ 
1.5 eV and fixed U = 3.5 eV is shown. Two main en- 
ergy regions can be identified: In the high energy region 
a; < — 4 eV the behavior of the self energy is governed by 
the correlation [/, which again can be explained by sum 
rules, cf. Ref. 46, On the other hand, between —4 eV and 
eV, the slope of changes decisively with increas- 

ing Ae. However, the latter part of the slope may again 
be explained by the arguments of Sec. |IVB[ since the 
values of limc^^o ^^{ioj) again depends approximately 
quadratically on the disorder Ae and agrees for small 
Ae. For Ae = 1.5 eV> t, the system is already insulat- 
ing, which is emphasized by an almost divergent in 
the vicinity of = 0. 
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FIG. 19: (color online) The imaginary part of the self energy 
^S(icj) on the Matsubara axis for T = 290 K, [/ = 3.5 eV 
and increasing values of the disorder potential As. (inset) The 
filling per spin nNa,nvac at the two sites CoNa,Covac, respec- 
tively, as a function of the disorder potential Ae. 



Another interesting numerical result is the filling of the 
diig orbital on the two sub lattices CoNa and Covac as 



a function of Ae, as visualized as an inset in Fig. 19 
The filling n^a rises in the analyzed Ae-interval basi- 
cally linearly to 1 as the disorder potential Ae increases, 
whereas nvac decreases. For very large disorder Ae = 1.5 
eV the fillings approach their maximum and minimal val- 
ues nNa=l and nvac=0-5, respectively. 
After the self energy, we investigate the spectrum for in- 
creasing disorder potential Ae. The k-integrated spectra 
A{uj) can be found in Fig. [20] (top). For Ae = ^ 1 eV 
the spectral weight at the Fermi edge decreases, but the 
asymmetry of the spectral distribution with respect to 
cj = increases. A further increase Ae > 1 eV leads to 
the formation of a band gap and a alloy band splitting 
created by disorder. 

Though less spectral weight implies decreasing ther- 
mopower contributions for both electrons and holes (see 
/Ctot in (28) and Figjl]), asymmetry in the spectrum with 
respect to the Fermi edge uj = enhances and the 
thermopower thus non-trivially depends on the disorder. 
In Fig. 20 (bottom), we show the numerical results for 
the thermopower 6* as a function of the disorder poten- 
tial Ae. The slope shows a flat maximum of the ther- 
mopower S = 80/iV/K for Ae ~ 0.7 eV. For larger 
Ae ^0.9 eV, transport gets more and more suppressed 
as the spectral weight is shifted away from cj = an al- 
loy band gap forms. Again, the k-resolved spectra and 
contributions to the thermopower shown in Fig. [23j en- 
hance the understanding of these effects: For no disorder 
Ae = eV, we obtain a weak renormalization as dis- 
cussed in Sec. IIVCI With a disorder of Ae = 0.7 eV 
all the charge carriers get scattered, which can be seen 
in a more diffuse k-resolved spectrum (corresponding to 
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FIG. 20: (color online) Top: k-integrated spectra for U = 3.5 
eV, T = 290 K and increasing disorder potential As. Spec- 
tral weight at the Fermi edge = diminishes with increas- 
ing Ae, but at the same time the asymmetry is increased. 
Bottom: The thermopower 5* as a function of the disorder 
potential Ae for a temperature T — 290 K. Starting from 
5* = 64 /jV/ K for Ae — eV the thermopower reaches a max- 
imum S = 79yiV/K for Ae = 0.9 eV. For larger Ae > 0.9 
eV, a continuous decrease of S can be observed as the system 
becomes insulating. 



a larger ^T^{iujm))' However, as the corresponding con- 
tribution to the thermopower reveals, electrons become 
scattered stronger than holes. Therefore the net enhance- 
ment for the transport property S in Fig. |20] (bottom) , 
can be explained by an increase of spectral weight for the 
holes UJ < compared to the purely disordered model, 
which is larger than the, also observed, increase of spec- 
tral weight for electrons {u > 0). But for Ae > 0.8 eV 
the contribution of holes to the thermopower starts to di- 
minish by a larger amount than can be accounted for by 
less negative contribution of the electrons. All charge car- 
riers then are scattered at a high rate, the thermopower 
S consequently drops as the system becomes insulating. 
After the numerical results for increasing disorder, we 
now fix Ae = 0.55 eV and gradually increase the cor- 
relation U = ^ 5 eV. In Fig. |2l] we plot the corre- 
sponding imaginary part ^T^{iuj). Surprisingly, the value 
lim^^o ^S(zco') first drops from 0.15 eV below 0.1 eV and 
then stabilizes around 0.1 eV for increasing correlation U. 



17 




iio (eV) 

FIG. 21: (color online) The imaginary part of the self energy 
^^{iuj) on the Matsubara axis for T = 290 K, Ae = 0.55 eV 
and increasing values of the correlation U. The black lines in 
the magnification are parabolic fits of the slopes. Note that the 
slopes converge to similar values for a; ^ 0, yet increase with 
larger U. For U < As, the value of lima;^.o ^^(icj) appears to 
be governed by the disorder, while for larger U the shape of 
^^(icj) is mainly determined by correlation effects. 



This indicates a rather non-trivial impact of a changing 
U on an already disorder system (the situation is further 
complicated since the Luttinger theorem does not hold 
due to disorder, and non-zero temperature, and the van- 
Hove like peak around the Fermi edge appears to react 
very sensitive to changes of U). 

The effect of increasing correlation on the k-integrated 

(top). A shift of spectral 
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spectrum is shown in Fig. 
weight from just above to just below the Fermi level 
can be observed. This indicates an increase of the ther- 
mopower for larger /7, which is indeed found as it is shown 
in Fig. 22 (bottom), where we show S function of U. 
The enhancement of 5 = 25 ^ 64/iV/K is extreme as 
U = ^ 3eV. However, the growth of S gets saturated 
for [/ > 5 eV. 

To further improve our understanding, we again investi- 
gate the k-resolved spectra and corresponding contribu- 
tions to the thermopower, cf. Fig. [24] Without correlation 
[/ = 0, there is little spectral weight directly around the 
Fermi edge. However, in that energy interval a larger U 
effectively packs spectral weight to larger extent for holes 
than for electrons. Thus, the drastic enhancement of S 
can be understood. For larger U ^ 5 eV, the change at 
the Fermi edge is not decisive, resulting in the satura- 
tion of the thermopower S with respect to U shown in 
Fig. [22];bottom). 



E. Effects of the doping 

The considerations of the previous sections are mostly 
thought experiments to identify the important effects in 
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FIG. 22: (color online) Top: k-integrated spectra for Ae = 
0.55 eV, T = 290 K and increasing correlation U. Spectral 
weight is shifted from right above to right below the Fermi 
level uj — 0. Bottom: The thermopower 5* as a function of the 
correlation U for a temperature T = 290 K. The thermopower 
increases from S ^ 25fiV/K for U — eV up to a maximal 
value of 5* = 67/jV/K for 1/^4 eV, where the effect it 
saturated. 



the model compound Nao.7Co02. On the other hand, 
the typical external parameters one can control in ex- 
periment, are the temperature T and the doping XNa- 
The gradual change or replacement of specific elements 
within a chemical compound has been used since ancient 
times. During the course of the 20th century, doping has 
become a major tool in technical applications and exper- 
imental materials science. 

In semiconductor applications, one usually dopes by re- 
placement: In addition to the main contents (e.g. Si), a 
similar element is added which leaves the crystal struc- 
ture in principle unchanged and takes the position of a 
fraction of the main chemical component (s). In the band 
structure, this procedure creates new bands which are 
energetically close to the ones of the desired charge car- 
rier, and can therefore be much easier thermally activated 
than in the pure compound. 

In contrast to replacement, Nao.7Co02 is doped by 
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FIG. 23: (color online) k-resolved spectra A{k,uj) (left) and corresponding thermopower spectral densities /Ctot (right) for 
various values for the disorder potential As = {OeV, 0.7eV, leV} at T = 290 K. The color code is given by the colorbar and the 
tight-binding fit is visualized as dashed line. The electronic contribution cj > to the thermopower is increased for As = 0.7 eV, 
but the hole contribution < is enhanced even more. For large values of As = 1 eV spectral weight is diminishing above and 
below the Fermi edge. 
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FIG. 24: (color online) k-resolved spectra A{k,(jj) (left) and corresponding thermopower spectral densities /Ctot (right) for 
various values for the correlation [/ = {0 eV, 2 eV, 5 eV} at T = 290 K. The color code is given by the colorbar and the 
tight-binding fit is shown as a dashed line. Spectral weight is shifted from above, directly below the Fermi edge cj = as 
U — ^ 2 eV, which increases the hole contributions to S while decreasing the electron contribution. This effect is even 
stronger for U = 5 eV. 
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FIG. 25: (color online) The thermopower S over the tempera- 
ture T for [/ = 3.5 eV, As = 0.55 eV and various doping XNa 
computed by DIMFT compared to experiment. Transport in 
xy-direction is assumed, if^y = Vx ~^ Vy. The thermopower 
increases non-linearly with temperature T as discussed in 
Sec. [IVAl For increasing doping XNa, the thermopower in- 
creases for all temperatures T. 



changing the sodium content, cf. Refs. [32H341 Since the 
single electron of the alkali Na becomes a valence elec- 
tron in the cobaltate C0O2 layers, this corresponds to 
shifting the Fermi energy Ep within the ai^ band of the 
non-interacting band structure (Fig. [6|) . 
As discussed in Sec. [ll] and in Ref. 11, the "pudding 
mold" -like slope of the ai^ band is expected to be a 
very important feature for a high thermopower S. The 
closer the Fermi energy Ep is to the flat top of the band, 
the larger the effect will be. Then, one charge carrier 
highly outmatches the other in terms of group velocity 
\/e{k). We therefore expect the thermopower to increase 
for higher XNa- In fact, as can be seen in Fig. 25, the com- 
putations yield a larger thermopower S as the doping XNa 
increases from 0.6 0.8. Here, we plot the thermopower 
S over temperature T for various doping x-^a and com- 
pare to experiment (see Ref. 9 for a corresponding com- 
parison of the resistivity). The experimental results are 
in good agreement with our data, even though a direct 
comparison is made difficult, by the (in this case) un- 
known orientation of the crystal in the experiment. The 
result of larger thermopower S with increasing can 
be qualitatively explained by band structure arguments. 
Still, there can be additional effects due to correlations 



and disorder. The k-integrated spectra in Fig. 26 indi- 
cate that spectral weight is shifted from electrons u > 
to holes < as XNa = 0.6 0.8. This further enhances 
the "pudding mold" effect of the group velocity Vs. The 



k-resolved spectra in Fig. 27 show another result: The 
"pudding mold" contribution of the k-path in the xy- 
plane T ^ K important for XNa = 0.6,0.7 is joined by 
the kz ^ contribution A ^ H diS the Fermi level moves 
down. That provides an explanation why the dampen- 
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FIG. 26: (color online) k-integrated spectra for U = 3.5 eV, 
As = 0.55 eV, T — 290 K and various doping XNa- For 
a^Na = 0.6 the Fermi liquid peak is the most pronounced as the 
filling of the site Covac is closest to 1/2. In the vicinity of the 
Fermi edge, spectral weight is shifted from cj > to cj < as 
XNa = 0.6 0.8. (bottom) The imaginary part of the self en- 
ergy ^ll{iuj) on the IMatsubara axis for U — 3.5 eV, As — 0.55 
eV, T = 290 K and various doping XNa = {0.6,0.7,0.8}. 
(inset) IMagnification around uj — and corresponding tan- 
gents (black). All slopes converge to similar values which are 
independent of the disorder XNa, but dependent on the disor- 
der potential As and U . 



ing of the linear increase of the thermopower S becomes 
smaller as XNa = 0.7 — ^ 0.8. 
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FIG. 27: (color online) k-resolved spectra A{k,(jj) (left) and corresponding thermopower spectral densities /Ctot (right) for 
[/ = 3.5 eV, As = 0.55 eV and various values of the doping XNa- The color code is given by the colorbar and the tight-binding 
fit is visualized by a dashed line. The (negative) electronic contribution to the thermopower is smaller as XNa = 0.6 0.8, 
thus the thermopower S shows a net increase. For XNa = 0.8 also the "pudding molds" between A ^ H and A ^ L becomes 
important. 
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V. CONCLUSION 

In Naa;Co02 two distinct effects enhance the hole cur- 
rent in comparison to the electron one: a larger group 
velocity and a larger number of holes or, more pre- 
cisely, a larger hole spectral weight^^, also see Fig. [l] 
As both effects point in the same direction, we obtain a 
large electron- hole imbalance. Consequently, a tempera- 
ture gradient results in a higher hole current, and, hence, 
a large positive thermopower S. The larger group veloc- 
ity of the holes can be understood already from a one- 
particle LDA point of view. It is caused by the particu- 
lar pudding-mold type of bandstructure^^. In this paper 
we analyzed the second mechanism which is a genuine 
correlation effect. The microscopic origin is a complex 
interplay of Coulomb repulsion U and a (disordered) po- 
tential Ae due to the Na vacancies. We have hence taken 
U and Ae as free parameters (not ab initio ones as in 
Ref. 9). The Coulomb repulsion alone (without disorder) 
has only a rather minute impact. This is due to the fill- 
ing of 1.7 electrons per Co site. With 0.3 holes in an 
otherwise filled Co ai^ orbital, chances to have 2 holes 
on the same Co site are very small. Hence the impact of 
the Coulomb repulsion, whose main effect is to suppress 
such configurations, is small. Without the vacancy po- 
tential, electronic correlations in Naa;Co02 would there- 
fore be very weak. Without Coulomb repulsion on the 
other hand, the vacancy disorder smears out the spec- 
trum and particularly enhances the resistivity less than 
the current-heat-current correlation function is dimin- 
ished, so that the thermopower actually decreases. The 
story becomes very different if Coulomb repulsion and 
disorder potential are combined. Then, the disorder po- 
tential gives rise to two inequivalent lattice sites. Since 
the number of holes in the aig exactly corresponds to 
the number of sites with a higher potential for the elec- 
trons, we end up with a situation close to two electrons 
on the Co sites with adjacent Na-ion and one electron 
on the Co site with adjacent vacancy. In this situation, 
electronic correlations can be very strong for the half- 
filled Co sites with an adjacent vacancy. Because of this, 
we get a strongly renormalized spectrum with a sharp 
peak at the Fermi level, see Fig. [9] The other sites, with 
adjacent Na-ions, on the other hand, have only a size- 
able spectral weight on the hole-side of the Fermi level. 
Taken together, we have more (quasi-) holes but nonethe- 



less sharp peaks in the spectrum. In this situation, the 
thermopower increases up to 80 jiV/K with interaction. 
For a larger value of the disorder potential, the Coulomb 
interaction can even (counter intuitively) change the be- 
havior from an insulator with a very small thermopower 
to a metal with a big one. 

For tracing down the energy and momentum of the 
quasiparticles responsible for positive and negative ther- 
mopower alike, we have analyzed the various contribu- 
tions to the thermopower separately. In Fig. [T2j we show 
for which k and uj we have a significant contribution from 
the group velocity, the spectral weight, and for which in- 
terval holes and electrons can be thermally activated. We 
also introduced a thermoelectric spectral weight which 
combines these terms and hence provides for a k- and uj- 
resolved analysis of the thermopower. The biggest contri- 
butions naturally stem from those /c-vectors which have 
the largest group velocity, albeit the quantitative contri- 
bution is changing in comparison to the non-interacting 
case. Most interestingly, also hole contributions appear 
which, without interaction, were not present at certain 
/c-points. Generally, electronic correlations change the en- 
ergies and momenta with big contributions to the ther- 
mopower quite substantially. 

As the vacancy (disorder) potential Ae is most crucial for 
the thermopower and the physical properties of NaajCo02 
in general, it is worthwhile to have a closer inspection 
of this parameter. Actually, in LDA calculations, Ae is 
rather small, i.e., ~ 0.05eV(Ref. 10). This value is way 
too small to provide for the aforementioned charge dis- 
proportionation and the resulting strong electronic cor- 
relatio ns on the Co sites with vacancy. Marianetti and 
Kotliail^^^EII, argued that LDA overestimates screening 
so that the actually Ae is much larger, i.e., ~ 0.55 eV. 
This is also the value we have considered in our paper as 
the physically relevant one. Besides the overestimation 
of screening, we think that the electron-phonon coupling 
is an important source for the effective increase of Ae 
by ~ /(^Q {g- electron phonon coupling strength; uq: 
phonon frequency). 
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